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Abstract 

We construct a recursion operator for the family of Narita-Itoh-Bogoyavlensky infinite 
lattice equations using its Lax presentation and present their mastersymmetries and bi- 
Hamiltonian structures. We show that this highly nonlocal recursion operator generates 
infinite many local symmetries. 



1 Introduction 

Bogoyavlensky constructed in |T] the following two families of integrable equations for arbitrary 
p£N known as an integrable discretisation for the Korteweg-de Vries equation: 



fc=l fe=i 

V V 

vt = v(Y[ v k~ J\v-k), (2) 

k=l k=l 

where both u and v are smooth functions of n £ Z and time t. Here we use the notations 

ut = dt(u), Uj = S 3 u{n, t) = u(n + j, t) 

and S is the shift operator. The same notation is also used for the function v(n, t) and later 
for the function w(n,t). Equation ([1]) is called the Narita-Itoh-Bogoyavlensky lattice due to the 
contributions of Narita and Itoh [2j [3] . 

Taking p = 1 in both ([1]) and ([2]), we get the well-known Volterra chain 

Ut = u(u\-U-\). (3) 

Thus these two families can be regarded as the generalisation of the Volterra chain. For fixed p, 
equation ([2D transforms into ([T]) under the transformation u = YVk=G v k ■ 

The so-called the modified Bogoyavlensky chain is given by 

p v 

w t = w2 (y\ w k - n w-k), (4) 

k=l k=l 

which is related to ([T]) by the Miura transformation u = Y\ I k=o w k- 
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All systems (JTj), ([2]) and can be defined for finite lattices, that is, the independent variable 
restricts to 1 < n < N for fixed N € N. Extensive work has been done to study their Hamiltonian 
structures, associations with classical Lie algebras, the r-matrix structure, etc. We refer the 
reader to the chapter on Bogoyavlensky Lattices in Suris' book |4J including the references 
mentioned in Bibliographical remarks in the end of this chapter. 

In this paper, we consider these systems defined on infinite lattices, i.e., n G Z and study 
their recursion operators and Hamiltonian structures for arbitrary p G N. We concentrate on 
system (pQ). The recursion operators and Hamiltonian operators for ([2]) and @ can be obtained 
via the Miura transformations mentioned before. 

A recursion operator for the well-known Volterra equation ([3]) is known to be 

» = uS + u + u 1 + uS~ 1 + ut(S -I) -1 - (5) 

u 

= u(l + S^Su - - S' 1 )' 1 ^ 1 . (6) 

It generates infinitely many local higher symmetries 
u tl = u(ui - 

u t2 = K(n tl ) = uu\(u + ui + u 2 ) - u-\u(u-2 + u-i + u) 



The Volterra equation ([3]) is a bi-Hamiltonian system 

Ut = H\5 u u = ^ 2 5 M ^, (7) 
where 5 U is the variational derivative with respect to dependent variable u, 

Mi = u{S-S' 1 )u, and U 2 = WHi = «(1 + - + S)u . 

When p = 2, system ([1]) becomes 

u t = u(u 2 + u\ — U-x - n_ 2 )- 

Its recursion operator, mastersymmetry and bi-Hamiltonian structure were presented in [5], 
where the authors also gave one Hamiltonian operator for arbitrary p G N as follows: 

H = u(£s i -Y J S- i )u. (8) 

i=l i=l 

Indeed, we can write system ([JJ as ut = H5 u u. 

Recently, Svinin [6] derived the explicit formulas of generalised symmetries for system (UJ in 
terms of a family of homogeneous difference polynomials by considering it as a reduction of 
a bi-infinite sequence of the KP hierarchy. The properties of these homogeneous difference 
polynomials [6] enable us to prove the locality of symmetries in Section HJ 

In this paper, we first construct a recursion operator of ([1]) for arbitrary p and then show that 
it indeed generates local symmetries. These symmetries are the same as the ones presented 
in [6] (modulo signs). Using the recursion operator, we further obtain the second Hamiltonian 
operator of ([T]), from which it naturally follows that the Narita-Itoh-Bogoyavlensky lattice is a 
bi-Hamiltonian system. 
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The arrangement of the paper is as follows: In section [21 we give the necessary definitions and 
fix the notation. The main result will appear in section O where we give the details on how 
to construct the recursion operator for all p G N via the Lax representation. In section [H we 
prove by induction that the recursion operator indeed generates local symmetries . We finish 
the paper with a discussion on recent results on the discrete Sawada-Kotera equation [7, 8j. 

2 Basic definitions of differential-difference equations 

In this section we introduce some basic concepts for differential-difference equations relevant to 
the contents of this paper. More details on the variational difference complex and Lie derivatives 
can be found in [TO] ■ 

Let u be a function of variable n G Z and time variable t. An evolutionary differential-difference 
equation of dependent variable u is of the form 



where [u] means that the smooth function K depends on u and its shifts. 

Assume that K[u] depends on a finite set of variables Uj, where q < j < p and q,p G Z with 
d Uq K ^ and d Up K ^ 0. We say K[u] is of order (q,p). The difference between p and q, p — q, 
is referred to as the total order of K. For example, the Volterra chain ([3D is of order (—1,1) 
with total order 2. 

All such smooth functions form a difference ring denoted A with the shift operator S as its 
automorphism. For any element / G A, we define an equivalent class (or a functional) //by 
saying that two elements f,g £ A are equivalent if / — g G Im(«S — 1). The space of functionals 
is denoted by A', which does not inherit a ring structure of A. 

A vector field (derivation) v is said to be evolutionary if it commutes with the shift operator 
S. Such vector filed v is completely determined by a smooth function P[u] G A. We call it the 
characteristic of vector field v. We present it as 



For any two evolutionary vector fields with characteristics P[u] and Q[u], we define a Lie bracket 
as follows 



where Dp = ^ ■ is the Frechet derivative of P. The evolutionary vector fields form a Lie 

algebra denoted by f). We simply say P G f). 

Definition 1. An evolutionary vector field with characteristic P[u] is a symmetry of system (0) 
if and only if [K, P] = 0. 

Equation ([9]) is said to be integrable if it possesses infinitely many higher order symmetries. 
Often the symmetries of integrable equations can be generated by recursion operators |llj . 
Roughly speaking, a recursion operator is a linear operator K : f) — > fj mapping a symmetry to 
a new symmetry. For evolutionary equation Q, it satisfies 



u t = K[u] 



(9) 




[P, Q] = D Q [P]-D P [Q], 



Dx[K] = [D K , 3?] , 



(10) 
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where -Dj{[-ftT] is the Frechet derivative of 9ft along the evolutionary vector field K. 

Symmetries of integrable systems can also be generated by mastersymmetries |12| . An evolu- 
tionary vector field r is a mastersymmetry of equation ([9]) if [r, Q] is a new symmetry whenever 
Q is a symmetry. Often the vector field r involves non-local terms. 

Recursion operators for nonlinear integrable equations are often Nijenhuis operators, that is, for 
any F 6 f) the operator 5? satisfies 

D U [SRP] - [Ar P , »] = 3ft (D n [P] - [Dp, »]) . (11) 

Therefore, if the Nijenhuis operator $ft is a recursion operator of ut = K, the operator is also 
a recursion operator for each of the evolutionary equations in the hierarchy ut = $l k K, where 
k = 0,1,2,... . 

Nijenhuis operators are closely related to Hamiltonian and symplectic operators. Their interre- 
lations were discovered by Gel'fand &: Dorfman |13t 114] and Fuchssteiner & Fokas |15| I16|. The 
general framework in the context of difference variational complex and Lie derivatives can be 
found in [91 [10] . A generalised definition of bi-Hamiltonian systems is given in |17l 118] . Here we 
recall the basic definitions related to Hamiltonian systems. 

We denote the space of functional (/-forms by Q q , q £N starting with 0° = A'. For any vertical 
1-form uj = h^ty on the ring A, there is a natural non-degenerate pairing with an element 
Pet): 



<u, p>= jY, hii)slp = j w2 s ~* h{i) \ p =<^«s^/i w , p>. 

i \ i / i 



(12) 



Thus the vertical 1-form uj is completely defined by £u = S l h^) u. We simply say ^ € SI 1 . 

The pairing (|12|) allows us to give the definition of (formal) adjoint operators to linear (pseudo-) 
difference operators [9l ITU], 

Definition 2. Given a linear operator A : t) — >• SI 1 , we call the operator A^ : t) — > SI 1 the adjoint 
operator of A if < AP, Q >=< A^Q, P >, where P, Q G h. 

Similarly, we can define the adjoint operator for an operator mapping from SI 1 to f), from f) to 
t) or from SI 1 to SI 1 . We say an operator A is anti-symmetric if A^ = —A. 

For any functional Jg E A', we define its difference variational derivative (Euler operator) 
denoted by 5 u (fg) € SI 1 as 



Definition 3. An anti- symmetric operator % : SI 1 — > f) is Hamiltonian if and only if the bracket 
defined on A' as follows: 



I 1, I 9 } = l {5uf mu9) 



satisfies the Jacobi identity. 

Definition 4. We say a Hamiltonian operator H : SI 1 — > f) is a Hamiltonian operator for 
equation (0|) if its Lie derivative along evolutionary vector field K 6 f) vanishes, that is 

L K H:= D H [K}- D K H-HD !f K = 0. (13) 
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Alternatively, we can define Hamiltonian operators in term of the Schouten bracket. 



Definition 5. Let TL,V : Q — > f) be two anti- symmetric operators. The Schouten bracket ofH 
and V is the trilinear mapping [[H, V]] : O, 1 x Q, 1 x {I 1 — > 0° defined by 



[[H, :=< 6,^[n 3 ](6) > + < iuD v [M?M2) > + < &,-D«P>£i](&) > 

+ < &,Drmi](&) > + < &,D H [VZ 2 KZi) > + < &,D v [H&](tii) > ■ 



By direct calculation, it can be shown that an antisymmetric operator % is Hamiltonian if and 
only if [[H, W\\ = 0, [HI 02]. Moreover, the relation between Nijenhuis operators and Schouten 
brackets can be formulated as follows: 

Theorem 1. Let H\^H.2 be two anti- symmetric operators and 3? = HxH^ ■ Then 



where P = %2£i and Q = Hi£,2- 

It follows that operator K = TiiTi^ 1 is Nijenhuis if %i and I-L2 form a Hamiltonian pair. 

Definition 6. An evolutionary equation (0|) is said to be bi- Hamiltonian if there exist two 
Hamiltonian operators TLq and Hi, and two Hamiltonians fo and f\ such that 



where 5 U is the difference variational derivative with respect to u. 

3 Construction of recursion operators from a Lax representation 

In general, it is not easy to construct a recursion operator for a given integrable equation although 
we have the explicit formula ()10p . The difficulty lies in how to determine the starting terms of 
3ft, i.e., the order of the operator, and how to construct its nonlocal terms. Many papers are 
devoted to this subject, see [2011271122] . 

If the Lax representation of an evolutionary equation is known, an amazingly simple approach 
to construct a recursion operator was proposed in [23] and later applied for lattice equations 
[24] . This idea can be used for Lax pairs that are invariant under the reduction groups [18] . 
In this section, we use the same idea to construct a recursion operator of system (pQ) from its 
Lax representation. The recursion operators of (|2|) and can then be derived via the Miura 
transformations. 

The Lax operator of ([T]) for any p E N is given in [I] as follows: 



< C,D^P](Q) - MDx[P](Q) - D^Q](P) + sR.Dsr[Q](P) > 

= \[[%x, Wi]](£,£i,6) + |[[H2, ft2]](» t2 £,£i,6)-[[fti, %]](Kte,6,6) 



(15) 



u t = Ho5 u f = %\8 u fi, 



L = S + uS~ p . 



(16) 



We can compute its hierarchy of symmetry flows using the formula 
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The idea to construct a recursion operator directly from a Lax representation is to relate two 
difference operators such as i?( ni ) and £>( n2 ) and to find the relation between the corresponding 
flows 1 231. 



We rewrite (|16|) into the following matrix form 

L(A) = S - AC/ (0) - C/ (1) :-- 



S-C(X), 



(17) 



where = (ufj) and = (u^j) are (p + 1) x (p + 1) matrices. The matrices and 
are given by 



/I 




o \ 




and C/ (1) 



••• 
\ • • • / 

respectively, that is, their non-zero entries are 

1 



1 










-u \ 




= 1, u\l i , = — ii and u„- 



1 
1 



1,2,--- ,p. 



11 — x > 1,P+1 ~~ " "i+l,* 
Further we rewrite the operator in a matrix form denoted by B(X). So the symmetry flows 
can be obtained by the zero curvature equation 



C(X) t = S(B(X)) C(X) - C(X)B(X). 



Make the Ansatz 



P+i 



B(X) = X P+1 B{X) + W(X), where W(X) = X p+1 -' l A^, 



(18) 



(19) 



i=0 



4?: 



are (p + 1) x (p + 1) matrices with the only non-zero entries being a^L 7 - for 



where .A® 

1 < i < P + 1- Here we read i + j as (i + j) mod (p + 1). For simplicity, we shall continue 
to denote the index I when I > p + 1 instead of / mod (p + 1). So both A( ) and A.( p+1 ) are 
diagonal matrices. 

The Ansatz W is invariant under the following transformation 

r : W(X) i — y PW(aX)p- 1 , 

where P is a diagonal (p+ 1) x (p + 1) matrix given by Pa = a 1 and a = e 2?rl// ( p+1 ) since we have 
P~ 1 A^P = ff'i". The transformation satisfies r p+1 = id and therefore generates the group 
Zp_|_i. The reduction groups of Lax pairs have been studied in [25| [26] . 

The zero curvature condition (j!8j) for B(X) leads to a formula for computing a recursion operator 
as follows: 



C tn+1 = S(B(X)) C(X) - C(X)B(X) = X p+1 C tn + S(W)C - CW. 



(20) 



Substituting (fT7|) and (fl~9j) into ([20]) and collecting the coefficients of this A-polynomial, we 
obtain 



A p + 2 : 
X p+l : 

^p+l—i 
A°: 



S(A^)UW -U^AW =0; (21) 
[/£) + 5(^( 1 ))L/(°) - + S(A^)U^ - uMA<® = 0; (22) 

S(A^)U^ - U {0) A^ + 5(AW)trW - C^WaW =0, 1 < i < p; (23) 



(24) 
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Lemma 1. Assume that p > 2 and 1 < i < p — 1. T/ie matrix equations 

S(A^)U^ - U^A^+V + S^C/W - E7WA» = (25) 

are equivalent to 

5(aS^J) + 5(a2 2i2 )-a« 4il = 0; (26) 

-°££-< + 5 ( a £ +2 - i ) + ^+i, P+ i- = °; ( 27 ) 

-«5(oiSi,i) -«S+i = 0; (28) 
«S(ag i+1)Z+1 )-ag M = 0, 2<l<p, l^p + l-i. (29) 

Proof. We directly compute the multiplications of matrices and write out their non-zero entries, 
respectively: 

(S(A^)U^) W = S(a£+>?? = 5(oS£}); 



,0 \„.M _ c/„W 



(5(A«)C/(D), +i+y = 5(a« +1 , W X^ = 5(aW +1) , +1 ) 5 1 < Z < p; 



) _ 



w, = «g> 4=4. i<i<p+i, i^p+i 



We are now ready to write out the entries for the matrix equations, which leads to the formulas 
stated in the lemma. ■ 

Specifically, when i = 0, using the proof of Lemma [T] we obtain the equivalent conditions for the 
matrix equation (j25j) as follows: 

S(4;j) + 5(a<g)-ag = 0; (30) 
-ag +1 + S(ag) + na^ lip+1 = 0; (31) 
5K% +1 )-4? =0, 2</<p. (32) 

In the similar way, the matrix equation (j25[) for i = p is equivalent to 

5(a^) - + S(ag>) + «ag lfl = 0; (33) 

-«5(ag Xil ) - ag +1 = 0; (34) 
5(agVi) " 4 P \l = 0, 2</<p. (35) 
Notice that formula pi]L ([HT ]) . (p ]> and ([33] ! are valid for p = 1. 

Using formula (|2"6"|) - ([2"9"j) in Lemma [TJ we can now find the relation between a^ 1 ] and af^ r 
Lemma 2. Assume that p > 2 and 1 < i < p — 1. W^e /iaue 

agj 1 } = 5- 1 (5 i u - «5 i_1, )- 1 (5 i « - n5 i - p - 1 )5(ag 1 J. (36) 
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Proof. First using formula (J29J), we can show that 



a« =5- 1 (aS2 +1 ). (37) 



Indeed, if we take I = p + 2 — i + rin (|29p . it follows that 

a l+r,p+2-i+r = $ ( a 2+r,p+3~i+r) ' < r < i — 2. 



Thus we recursively obtain (j37|) . that is, 

,(*) _ e/^W \ _ _ ci-l /■„(*) 



l,p+2-i - °V a 2,p+3-J - • • • - O l a i,p+lJ- 

From ([28]) . it leads to a^ +1 = — u^a^ x ). Substituting it into to (|37p. we get 



(0 — c*-i», c 

l,p+2-i ■ 



5 l -^«S(a^ M ). (38) 



Similar as the proof of (|37|) . by taking / = 2 + r, where < r < p — 2 — i in (|29p when p > 2, we 
recursively show that 



a 



2+i,2 



SKy—^ 1 -^). (39) 



Now we substitute (f3"9"j) into (f2l)j) and it leads to 

<%£3) + ^("p+w-J - a il,i = °- ( 4 °) 

Formula ([38]) is valid for different value of 1 < i < p— 1. Thus we have a^^^j = — 5 l u5(a^ 1 ^). 
We substitute it and (f3"8"j) into (f2"T|) and it becomes 

^(ogSji) " 5*«5(«S?i fl ) + ua pti,p+i- = o- ( 41 ) 

From (j4"0j) and (JHJ), we eliminate a^|_ lp+1 _ i and obtain the relation between a^ 1 ] and 1; 
that is, 

af+V = (S*uS - uS^rHs'uS - uS^Ka^), 

which is equivalent to (j36|) as written in this lemma. ■ 



From (|34J) and (|35j) . we can see that (|38j) is also valid for i = p. Substituting it into (|33j) . we 
obtain 

(S - l)a^ +1 > = -5(ag) - iiag^ = (5^5 - u)^). 

Therefore, 

aft^ = (5 - l)- 1 ^ - u«S- 1 )5(ag 1)1 ). (42) 

We are now ready to compute the recursion operators for any p£N. 
Theorem 2. A recursion operator of Narita-Itoh-Bogoyavlensky lattice (CJ) is 

» = tt(5 - 5- p )(5 - I)" 1 ]J(5 P+1 -^ - uS^S^u - uS^y 1 ■ (43) 

8=1 



Since the difference operators are not commuting, here we use the notation f[i=i t° denote the 
order of the value i, from 1 to p, that is, Ili=i a * = a i a 2 ■ • • 

Proof. First from (|32p . we can show that 

og) = S^-\a£l ljP+1 ), 2</<p + l. (44) 

The next identity ([21} leads to (5-1) — 0. Here we choose the solution — 0, which 
makes it possible to find the relation between ut n+1 and ut n ■ We now substitute this into ([30]) 
and using ([44"]) we obtain 

og = -ag = -5^(4°i, P+1 )- (45) 
It follows from ((22J) and ((3TJ) that 

(1) , (0) n 

-«*n - °l,p+l + M «p+l,P+l = - 

We now show that a^p+i = —uS(a^\) for all p G N. When p > 2, we have o^p+i = ~~ ^(^l) 
by taking z = 1 in (|28|) . With p = 1, this is a result from ()34|) . Therefore for all p G N, using 
([43]) we have o^p+i = u< ^ p ( a p+i p+i) • This leads to 

(0) 



+ u(S* - l)a^ lip+1 = 



Therefore, we have 



Form (ED, we find that 



a^v^a-^v- (46) 



-u tn+1 = -«5(o^ ; ) + ua; +lp+1 , 
>+i) a _ „(p+i) 



o{a l+l l+1 ) - a u , l s « S P • 



This implies that 



Using (|42|) we rewrite it as 

u tn+1 = u(S - S->)(S - lyHS^u - n5~ 1 )5(ag lil ). 
Using Lemma [2] when p > 2, we obtain 

-Kp-1) 

Uin+i =n(5-5- p )(5-l)- 1 (5 p -u-n5- 1 ) [] (^m-mO^G^u-uS-*- 1 ).?^!) . 

i=i 

This is also valid for p = 1 under the convention that the empty product is equal to 1. Finally, 
using (|45|) and (|46p . we get the relation between ut n+1 and which gives rise to the recursion 
operator as stated in this theorem. ■ 

By directly checking the identity ([lip , we can prove that the operator K in Theorem [2] is 
Nijenhuis. Here we omit the proof since we have not found a neat way to present it in reasonable 
length even for a given p, 

Notice that equation ([I]) admits a scaling symmetry u, that is, [u, ut] = Uf. Knowing the 
recursion operator, one can obtain its master symmetry [12] . 
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Corollary 1. A master symmetry for equation is r = $tu, where 3? is the recursion operator 
given by Theorem^ 

Example 1. For p = 1 in ([7J it becomes the Volterra chain (0). From Theorem^ we get the 
same recursion operator as given by |5|). It follows from Corollary [I\ its mastersymmetry is 

t = K p= iu = nu{u\ — U-i) + u{u\ + u + 2u_i). 

The Lie derivative of its Hamiltonian operator Hi defined by ^7ty along vector field r leads to 

L T Hi = D Hl [r] - D T Hi - H X D\ = -H 2 , 

which is defined by It is also a Hamiltonian operator of the Volterra chain. 

For equation (H|), one Hamiltonian operator H given by (|8|) is known [5]. We now compute 
H = Wi, that is, 



u 



n = u(J2 s ~i 1 [ {sf+^u - us-^s^u - us-*)- 1 {Su - us- p ) si 

i=0 \ i=l ) \i=0 J 

This operator is anti-symmetric since it satisfies H) = —H. 
Proposition 1. The above anti- symmetric operator H is Hamiltonian. 

Proof. Here we sketch the proof. First by direct calculation we have L T % = —%. Since % is 
Hamiltonian, its Schouten bracket vanishes, that is, [[H, H]] = 0. It implies that [[L T H, H]] = 0, 
[2], implying [[H, H]] = 0. Due to K being Nijenhuis, it follows from ([15]) the Schouten bracket 
of H vanishes. Thus it is a Hamiltonian operator. ■ 

Hence the Narita-Itoh-Bogoyavlensky lattice is a bi-Hamiltonian system. Indeed, we can write 
equation ([1]) as 

In ii 

iH = H5 u u = mi8 u . (47) 

p+ 1 

The recursion operators and Hamiltonian operators for equations ([2]) and Q can be obtained 
by the Miura transformations. Notice that 



/p—l \ /p— l 



D u = u 



yj-s k ) = u ys k )- = u y—s k ) = u (ys k \-. 



Using Theorem [2j we obtain the following results: 
Corollary 2. For any p € N, a recursion operator of equation fcjQ) is 

— >p— l 

R v = v (S p - 1)-\S - S~ p ) ■ fl (<S P+1 ~^ - uS'^S^u - uS- 

i=i 

■{Su-uS- p )S p {vS -v)-\ (48) 
where u = 131=0 v k ■ A recursion operator of equation ^ is 

U w = wS- p - Y[ (S p+1 ~ l u - uS-'XSf-'u - uS- 1 )- 1 
i=i 

■(Su - uS- p )(l - S-f)- 1 ^ 1 - l)(wS - w)-\ (49) 

where u = IIaUo w k> 



10 



The Hamiltonian operators for equations ([2]) and Q are given in [5j. They are 




) 



) 



v 



and 



D- l %Di = wi- 
w{S - 1){S P - 



1-S- 1 5-1 



■w 



S p+1 - 1)~V 



w 



respectively. 



The equations ^ and @j are also bi-Hamiltonian systems. They can be written as 





w t = HwK || Wh = fcwHviSw In w 

k=0 



respectively. 

Here we make two remarks on the recursion operator 3? in Theorem [2) 

1. Every pseudo-difference operator can be uniquely represented by its formal Laurent series 
|27j . The highest order of the Laurent formal series of recursion operator 3? in Theorem [2] 
is uS p for all p£N. Thus the fractional power of this pseudo-difference operator does not 
exist with coefficients in u and its shifts, that is, we can not formally find an operator Q 
such that Q p = R. 

2. From the difference analog of the Adler Theorem [27} I28j. it follows that the residue and 
the logarithmic residue I > are canonical conserved densities of the Narita-Itoh- 
Bogoyavlensky lattice. 

4 Locality of symmetries 

Recursion operators often have nonlocal terms. One important question is whether the operator 
is guaranteed to generate local symmetries starting from a proper seed. Sufficient conditions 
for weakly nonlocal [29] Nijenhuis differential operators are formulated in |21l 130] : which are 
also valid for weakly nonlocal Nijenhuis difference operators [10]. This result is generalised 
to Nijenhuis operators, which are the product of weakly nonlocal Hamiltonian and symplectic 
operators [18] . However, the recursion operator given in Theorem [2] are not weakly nonlocal. In 
this section, we are going to directly prove the locality of symmetries by induction. 

To do so, we first introduce a family of homogeneous difference polynomials of degree I with 
respect to the dependent variable u and its shifts 




(50) 
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where k > 0, 1 > 1 and p > 1 are all integers. In particular, for any p € N we have 

p(i,k) = ^ Uj and P (z ' 0) = uu p u 2p ■ ■ ■ u {l _ l)p . (51) 

j=0 

Equation ([1]) can be written in terms of these polynomials as follows: 

ut = u{S - S~ P )V {1 ^ 1] . (52) 

This family of polynomials was first defined in [6], where the author amazingly gave the explicit 
expressions for the hierarchy of symmetries of equation ([1]) in terms of it. We are going to show 
that the recursion operator ()43[) generates the same symmetries as in this paper modulo a sign. 

It is easy to see that the polynomials (f5Q|) possess the following properties [6]: 

v (l,k) _ r (l,k-i) = UkS P{V {l ~ l ' k) ); (53) 
-pW _ s^fc-i)) = U(J _ 1)p 7?(J-i.*). (54) 

These immediately lead to 

( 5 - = Uk+l S p+1 (V«- 1 ^) - u {l _ 1)p V^ k \ (55) 

We now prove another important property. 
Proposition 2. For all l,p G N, we have 

(S^u - uS-^S-^V^+Vp-^ = (S^u - u s-^)S- lp+i+1 V il ' {l+1)p ~ i ~ 1) , < i < p. (56) 

Proof. Let us compute the difference between the left-hand side and the right-hand side of the 
identity ([56]) using the properties ([53]) and (Ell) : 



u p -iS~ lp+p (pV'V+Vp-V - S p(i,(i+i)p-i-i)^ _ uS -Ip ^p(l,(l+l)p-i) _ p(i,(l+i)p-i-i) 

u n .,S~ lp+p (u^ 1)p V^ l+1 ^) - uS~ lp ( U{l+1)p _ i spr^ l+1 ^) 



*-p— «<■ 

= u p . iU S- lp+p V (l -^ l+1)p - l) - uup^S-^V^ 1 ^ 1 ^ = 0. 

We proved the statement. ■ 
Notice that we can rewrite the recursion operator ()43|) in the form 

5ft = u(S - S' P ){S - l)- 1 {S p u - uS' 1 ) ■ \\ {S p ~ l u - uS-^iS^u - n5~ (i+1) ) 

i=i 

■{u-uS- p )- 1 . (57) 
Using it, we are able to prove the following result: 

Theorem 3. For Narita-Itoh-Bogoyavlensky lattice ([7]), starting from the equation itself, its 
symmetries Q l = ^R l (ut) generated by recursion operator |^3| ) are local and 

Q l = $( Ut ) = n (l - S-(P+l))sl-fr7>('+l.(l+l)P-l) for all < I 6 Z. (58) 
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Proof. First the statement is clearly true for I = from (|52p . Assume the statement is true for 
I — 1 > 0. Let us compute the next symmetry Q l . Taking i = p in (|56|) and using the induction 
assumption, we have 

(u - uS- p )S' lp+p V^ lp) = (u- uS-( p+l) )S- lp+p+1 V {l ' lp - 1] = Q 1 ' 1 . 

Hence 

Q l = ^(Q 1 - 1 ) = u(S - S- p )(S - ly^SPu - uS- 1 ) 

• Y[ (S^u - uS-'TH^u - u<S~ (m) ) ■ S~ lp+p V {hlp) , 

8=1 

where we used formula (|57p . We now recursively apply formula (|56p for i from p — 1 to 1 and 
obtain 

Q l = ^(Q 1 ' 1 ) = U (S - S~ P )(S - l)" 1 ^ - uS- 1 )S- lp+1 V^ lp+p - 1) 
= u(S - S- p )(S - I)- 1 (u p S- lp+p+1 V {l ' lp+p - 1) - uS- lp V {l ' lp+p - l) ^ 

= U (i - s-^s^is - i)- 1 (u {l+1)p s p+1 v^ l+1 ^ - u lp r^ l+ V p -v) 

= u(l - S~ ( - p+1) )S 1 ~ lp 'P {l+1 ' {l+1)p - 1) . 
Here we used formula ()55p for I being / + 1 and k being (7 + l)p — 1, that is, 

(5 - i)pP+i.0+i)p-i) = U{l+1)p S p+1 (V^ l+1 ^) - u lp V^ l+ ^ p - l \ 
We completed the induction proof of the statement. ■ 

5 Discussion and further work 

Recently, a family of integrable lattice hierarchies associated with fractional Lax operators was 
introduce by Adler and Postnikov [TJ [8] . One simple example is 

ut = u 2 (u p ■ ■ ■ u± — u_i • • • U-p) — u(u p -i ■ ■ ■ u± — n_i • • • ui- p )), 2 < p G N, (59) 

which is an integrable discretisation for the Sawada-Kotera equation. Notice that equation (|59p 
is a combination of equations (j3|) and ([2]) with different p. It can be considered as inhomoge- 
neous generalisation of the Bogoyavlensky type lattices. To study its algebraic and geometric 
structures, we need first know the answers to the Bogoyavlensky type lattices. This is the main 
motivation of this paper. 

As pointed by the authors of [8], an important open question of this new family of integrable 
equations is to construct their Hamiltonian structures. The present paper does not offer an 
immediate solution and further research is required. 
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